Optimal Entanglement Enhancement for Mixed States 
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We consider the actions of protocols involving local quantum operations and classical communi- 
cation (LQCC) on a single system consisting of two separated qubits. We give a complete description 
of the orbits of the space of states under LQCC and characterise the representatives with maximal 
entanglement of formation. We thus obtain a LQCC entanglement concentration protocol for a 
single given state (pure or mixed) of two qubits which is optimal in the sense that the protocol 
produces, with non-zero probability, a state of maximal possible entanglement of formation. This 
defines a new entanglement measure, the maximum extractable entanglement. 
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I. INTRODUCTION 



Entanglement is a basic quantum communication re- 
source which can usefully be manipulated to suit par- 
ticular tasks |l]j|]. In this paper we investigate the ma- 
nipulation of a single entangled mixed states comprising 
two separated single qubit subsystems. We consider two 
parties, Alice and Bob, who each control one subsystem, 
and who are restricted to carrying out local quantum 
operations and classical communication (LQCC). Specif- 
ically the quantum operations Alice and Bob are allowed 
to perform are local unitary transformations and local 
filtrations. The restriction to local quantum operations 
ensures that entanglement is indeed treated as a resource: 
if non-local quantum operations were allowed, Alice and 
Bob could create entanglement between them from ini- 
tially non-entangled states. 

The interest of this problem is that since any real- world 
quantum communication channel will be imperfect, even 
if Alice could create perfect maximally entangled states, 
she would never be able to share such states with Bob 
simply by sending one subsystem through the channel. 
So it is natural to ask whether Alice and Bob can use 
LQCC to obtain states with better entanglement from 
imperfectly entangled states. 

Various entanglement purification protocols have been 
suggested. If Alice and Bob share a number of copies of 
an imperfectly entangled known pure state, they can ob- 
tain a number of maximally entangled states by carrying 
out operations on each state individually or by collective 
operations on a number of shared states m. The col- 
lective algorithm has a higher asymptotic yield of maxi- 
mally entangled states in the limit in which the number 
of shared states tends to infinity. Efhcient collective al- 
gorithms which give a non-zero asymptotic yield of max- 
imally entangled states from entangled mixed states have 
also been described . 



In practice, though, the number of states will always 
be finite, and Alice and Bob will effectively share a single 
entangled state of two subsystems whose state spaces are 
finite dimensional. For this and other reasons — for ex- 
ample, Alice and Bob might actually have only one copy 
of an entangled state of some simple system, or it may 
be technologically difficult to implement collective oper- 
ations — it is interesting to see what Alice and Bob can 
achieve by gambling with the entanglement of a single 
state. That is, we would like to know how far the entan- 
glement of single states could be increased by LQCC if 
the outcomes of Alice and Bob's local measurements were 
favourable. The Procrustean algorithm of |l|] provides an 
answer to this question in the case of pure states. Here 
we answer the question for two qubit mixed states, and in 
the process illustrate a general approach to the problem 
based on identifying quantities invariant under LQCC. 

Though most mixed state entanglement distillation 
protocols discussed so far involve collective operations 
on many states, it has been established that there exist 
entangled mixed states for which single-state LQCC pro- 
tocols can increase entanglement [g[||. Conversely, it is 
known that there exist entangled mixed states, including 
the important case of the Werner states, for which no 
single-state LQCC protocol can increase entanglement 

We give here a complete description of the effect of 
LQCC on entanglement of a single copy of an arbitrary 
mixed state p of two qubits. It has been shown previously 
that if Alice and Bob's local density matrices are com- 
pletely random, they cannot increase the Entanglement 
of Formation (EOF) by LQCC ||-||. Here we show that 
if the local density matrices are not random and if the 
EOF is non- vanishing, Alice and Bob can always increase 
the EOF. Moreover we construct a procedure that max- 
imises the EOF of the final state. This procedure, which 
is unique up to local unitary transformations, leaves Alice 
and Bob with completely random local density matrices. 



II. MAIN RESULTS 

Throughout, the states considered are those of a single 
system comprising two separate single qubit subsystems. 
We use the following facts. [g[ 

1. The LQCC protocols we consider map the state p 
to states of the form 



Tr(A®Bpy4t (giSt) ' 



(1) 



where A and B are arbitrary operators that act on 
Alice and Bob's Hilbert space respectively. The 
only condition they must obey is ^^A < /2, 
B'I'B < l2- The protocol succeeds with probability 
lT{A®BpA'^ ®B^). We need not consider the most 
general local protocols in which the final state con- 
sists of mixtures of states of the form eq. (0) since 
mixing decreases the EOF. 

The operators A and B can be written as 



A®B^ UAf^""' ® UbI' 



(2) 



where Ua^ Ub are unitary and the filtrations / are 
defined by 

(3) 

Here p, v, a, b are real numbers, /„ denotes the iden- 
tity operator in n dimensions, and the vector cr — 
{cTi, (J2, era} has the Pauli matrices as components. 
We can also write these operators as ^ = UaFaU'^, 

ai ' 



where Fa takes the form _ with the a; 

V "2/ 

real, < a,; < 1 and Ua, U'j^ unitary; similarly 
B — UsFbUq. We can thus write any non-trivial 
LQCC (i.e. any LQCC which is not the zero map) 
in the form 



iUa 



1 

UA 



U'a®Ub 



1 
Ub 
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B J 



(4) 



where 7 is a scale factor in the range < 7 < 1 
and < a A, olb < 1- 

2. The entanglement of formation (or EOF) of a pure 
state IV') is defined as -E(V') = —TrpA^^PA — 
-TrpBliapB where pA = TrB\4'){'>P\, Pb = 
TrA\4'){'4'\ S'l'e the local density matrices seen by 
Alice and Bob. For a mixed state the EOF is de- 
fined as H]: E{p) = min J2iPiF{'>pi) where the min- 
imum is taken over all decompositions of p into pure 
states p = Y^iPt\ipi){'ipi\- 

In the case of a mixed state comprised of two sin- 
gle qubit subsystems, Wootters [p3 has given an 



explicit formula for E{p), verifying an earlier con- 
jecture of Hill and Wootters Q. Let p = (J2 ® 
<J2P*<J2 ® (^2- Call Ai the positive square roots of 
the eigenvalues of the matrix pp written in decreas- 
ing order. Define the concurrence by 



C{p) ~ max(0, Ai — A2 — A3 — A4) . 



Then the EOF of p is 



E{p) = H{ 



l + ^l-C^ip). 



(5) 



(6) 



where H{p) = -plogap- (1 - p) log2(l -_p). 

3. Consider a general density matrix p of two qubits. 
It can be written as 

p = -{h + a.cr ®l2 + h® P-o- + Rijdi ® aj) . 



(7) 



In [|1 it was shown that under LQCC of the form 
eq. (0) the positive square roots of the eigenvalues 
of the matrix pp transform as 



Ai — > Aj — — — : ^Ai 



(8) 



t{p;p,a,m;i^,b,n) 

where t is the probability that the LQCC succeeded 

t{p;p,a,iLn;i^,b,n) ^ 

p^iy^ [(I + a2)(I + b'^) + 2a(I -I- 6^)11 ■ a + 
2b{l + a'^)in- iS + AabRijUimj]. (9) 

Thus the concurrence also transforms as 



t{p;p,a,m;v,b,n) 



(10) 



It follows from eq. (||) that the ratios \i/Xj are 
invariant under LQCC. We add here the necessary 
qualification that the LQCC must be invertible. 

Now our argument runs as follows. We consider states 
p which have non-zero EOF and which are not Bell diag- 
onal (recall that a state is Bell diagonal if all its eigen- 
vectors are maximally entangled; equivalently it satisfies 
trA{p) = trB{p) — ^l2, ie. a = (3 — in the expres- 
sion for p given in eq. (H)). We show in Theorem 1 that 
there is an LQCC protocol which increases the EOF of p 
with non-zero probability. We show further in Theorem 
3 that this process can be iterated to obtain an LQCC 
protocol which, with non-zero probability, maps p to a 
Bell diagonal state with maximal EOF. In Theorem 4 we 



show that this is the unique optimal protocol up to local 
unitary rotations. 

Theorem 1 Let p be a density matrix of a state 

with non-zero EOF written as in eq. (R). If a or /3 are 
non-zero, then there is an invertible LQCC A (^ B map- 
ping p with non-zero probability to a density matrix p' 
with higher EOF than p. 
Proof For small a and b, eq. (|l^) takes the form 



Cip') 



1 



1 + 2am. a + 26n/3 



C(p). 



(11) 



then there exists an invertible LQCC which, with non- 
zero probability, maps p to a Bell diagonal density matrix 
p' which has the maximum EOF of any density matrix 
obtainable from p by LQCC. 

Proof Since by Lemma 2 the space of normalised 

LQCC which do not decrease the EOF of p is compact, 
and the EOF is a continuous function, the lowest up- 
per bound on the attainable EOF is attained by some 
LQCC. The corresponding density matrix p' must have 
a' = /3' ~ 0, otherwise, by Theorem 1, its EOF could be 
increased. 



Hence if a or /3 are non-zero and if C{p) is non-zero we 
can always find an LQCC which, with non-zero probabil- 
ity, increases the EOF, by choosing appropriately small 
a and b and suitable m and n. 

We now need a technical lemma about the topology of 
the space R of LQCC operations which do not decrease 
the EOF of a given p. The result, namely that R is com- 
pact, is needed in Theorem 3. 

Lemma 2 Let p have non-zero EOF. There exists a 

positive bound 6{p) such that if the state p' has greater 
EOF than p and can be obtained from p with non-zero 
probability by LQCC, then there exists some LQCC from 
which p' can be obtained from p with probability greater 
than 5{p). Furthermore let R be the space of LQCC 
which succeed with non-zero probability in producing a 
density matrix with EOF greater than or equal to that 
of p. Then R is compact. 

Proof Fix p. If we write A ^ B in the form (0), p' 

is independent of the scale factor 7, so that any p' ob- 
tainable from p can be obtained by a normalised LQCC, 
taking the form (0) with 7 = 1. For e > 0, define S^ 
to be the set of normalised LQCC of the form (0) with 
min{aA, cts} — £• Let E'^ be the maximum EOF of any 
density matrix p' obtained from p by the action (Q) for 
some A(Ei B in S^. Since E'e is continuous in e and tends 
to zero as e tends to zero, there is some positive ep such 
that E^ is less than or equal to the EOF of p for e < eg 
and such that eo(p) is maximal with this property. Let 
T^g be the union for 1 > e > eo of Se- Now if a non-trivial 
LQCC A'S)B annihilates p, i.e. A^BpA'' <S)B^ = 0, then 
A® B\ipi) = for all i (where jV'i) are the eigenvectors 
of p with non-zero eigenvalue) . Thus either A or B must 
be a rank one projector up to a scale factor. Hence no 
LQCC in TgQ can annihilate p. Also T^g is compact. So 
the probability Ti(A (g) BpA^ (g) i?^) of obtaining p' from p 
via the LQCC is non-zero everywhere in T^^ and attains 
a non-zero lower bound 5{p) on the set. This is a lower 
bound for all LQCC increasing the EOF of p, since no 
LQCC outside T^g does. The compactness of R follows 
since it is a closed subset of T^g . 

Theorem 3 Let p written as in eq. (M) be a den- 

sity matrix with non-zero EOF. If a or /3 are non-zero, 



Theorem 4 Let p be the density matrix of a state 

with non-zero EOF. Then the Bell diagonal state p' which 
can be obtained from p by LQCC is unique up to local 
unitary transformations. This p' has maximal possible 
EOF. 

Proof We start by calculating the positive square 

roots \i of the eigenvalues of the matrix pp. We order 
them as Ai > A2 > A3 > A4. The ratios -^ are invariant 
under the actions of invertible LQCC, see eq. (g). We 
characterise these ratios by the three numbers Ci — A^/Ai, 
i = 2,3,4. 

From Theorem 3 we know that p can be brought to Bell 
diagonal form by LQCC. We shall now show that the Bell 
diagonal form is uniquely specified, up to local unitary 
transformations, by the ratios Ci. To this end consider a 
Bell diagonal state pn = \{Ia-\- Rijcn <Jj) with positive 
EOF. Local unitary operations Ua €5 Ub transform pR to 
PR' = j{h + R[j(T, ®aj), where R' = (Oi)^i?(02) for 
some elements Oi and O2 of S0(3): any pair of Oi can 
be produced by suitable choices of Ua,Ub- By using a 
singular value decomposition |ll[| of R, we can find or- 
thogonal Oi such that R' is diagonal, so we can find local 
unitary operations mapping pn to the form 



1 



3 



Pri,r2,r3 = T ( A + '^nai ® (Ti) 

4=1 



(12) 



with all ri having the same sign and with ri < r2 < r^. 

Now Pri.r2,r3 = Pri,r2.r3, hcncc the eigenvalues of 
Pri.r2,r3 are equal to the A^. These eigenvalues are 
jil-ri-r2- rg), i(l + n + r2 ~ rg), i(l + r2 + r^ - 
ri), -1(1 + r^ + ri — r2). Since Pri,r2,r3 is assumed to 
be entangled, the r^ are all less then or equal to zero. 
(This may be verified by checking that when the ri are 
all positive the concurrence vanishes). We can now ex- 
press the ratio's Ci in terms of the r^. For instance 
C2 = (1 + ^2 4- r3-ri)/(l-ri -r2-r3). It is straightfor- 
ward to verify that the r^ can be uniquely expressed in 
terms of the Ci by inverting these equations. Therefore 
the Bell diagonal state of the form eq. (|l^) to which p can 
be brought is unique. 



III. CONCLUSIONS 



We have shown that any entangled state p of two qubits 
whose local density matrices are not completely random 
can be brought by LQCC to a unique (up to local unitary 
transformations) Bell diagonal state. No other LQCC 
can bring p to a state with more entanglement. To ob- 
tain an explicit expression for this optimal protocol, one 
should write explicitly the conditions that the density 
matrix p' obtained from has completely random local 
density matrices Ttap' — Trsp' = h- We have shown 
that these equations have a unique solution for the co- 
efficients a, m, b, n of the filtrations /^'°'™, f'^,b,n jj^ gqg^ 

(li- 

Our optimal protocol should be compared to the Pro- 
crustean algorithm for concentrating pure state entangle- 
ment of m which brings a non maximally entangled pure 
state to a maximally entangled pure state by LQCC. The 
main difference between the two methods is that the op- 
timal mixed state protocol generally requires Alice and 
Bob to carry out different filtrations and then tell each 
other whether the filtrations have succeeded. Only if 
both succeed do they obtain (and know that they have) 
the state with maximum cxtractable entanglement. The 
Procrustean algorithm on the other hand can be realised 
without classical communication between Alice and Bob, 
or with only Alice carrying out the filtration and com- 
municating the result to Bob. 

In H it was noted that the ratios Ci = A^/Ai arc in- 
variant under invertible LQCC. The argument used in 
proving Theorem 4 also shows that for entangled states 
they consitute an exhaustive set. Indeed we can bring 
any entangled p to the form eq. (O) which is charac- 
terised by three parameters r^ and they are in one to one 
correspondence with the Ci . This gives a characterisation 
of locally equivalent entangled density matrices. 



Our method also introduces an interesting combination 
of these invariants: the maximal extractable entangle- 
ment of a density matrix. This quantity provides a new 
characterisation of the entanglement of a state. It has 
the important property that it decreases under mixing 
(this follows from the convexity of the EOF [|l). 
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